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ABSTRACT 

We  present  results  of  a  high  angular  resolution  survey  of  massive  OB  stars  in  the  Cygnus  OB2  association  that 
we  conducted  with  the  fine  guidance  sensor  IR  (EGSlr)  on  the  Hubble  Space  Telescope.  EGSlr  is  able  to  resolve 
binary  systems  with  a  magnitude  difference  AV  <  4  down  to  separations  as  small  as  (y.'Ol.  The  sample  includes 
58  of  the  brighter  members  of  Cyg  OB2,  one  of  the  closest  examples  of  an  environment  containing  a  large  number 
of  very  young  and  massive  stars.  We  resolved  binary  companions  for  12  targets  and  confirmed  the  triple  nature  of 
one  other  target,  and  we  offer  evidence  of  marginally  resolved  companions  for  two  additional  stars.  We  confirm 
the  binary  nature  of  1 1  of  these  systems  from  complementary  adaptive  optics  imaging  observations.  The  overall 
binary  frequency  in  our  study  is  22%  to  26%  corresponding  to  orbital  periods  ranging  from  20  to  20,000  yr. 
When  combined  with  the  known  short-period  spectroscopic  binaries,  the  results  support  the  hypothesis  that  the 
binary  fraction  among  massive  stars  is  >60%.  One  of  the  new  discoveries  is  a  companion  to  the  hypergiant  star 
MT  304  =  Cyg  OB2-12,  and  future  measurements  of  orbital  motion  should  provide  mass  estimates  for  this  very 
luminous  star. 

Key  words:  binaries:  visual  -  open  clusters  and  associations:  individual  (Cyg  OB2)  -  stars:  early-type  -  stars: 
massive  -  techniques:  high  angular  resolution 
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1.  INTRODUCTION 

Massive  stars  (>10Mq)  play  a  fundamental  role  in  the  evo¬ 
lution  of  the  universe,  from  influencing  galactic  dynamics  and 
structure  to  triggering  star  formation  through  their  spectacularly 
violent  deaths.  It  has  been  well  established  that  massive  stars 
have  a  higher  binary  frequency  than  lower  mass  stars  (Mason 
et  al.  2009;  Raghavan  et  al.  2010),  and  the  binary  fraction  may 
be  as  high  as  100%  for  massive  stars  in  clusters  and  associations 
(Mason  et  al.  1998,  2009;  Kouwenhoven  et  al.  2007;  Chini  et  al. 
2012).  Because  most  massive  stars  are  born  in  clusters  and  asso¬ 
ciations,  their  multiplicity  properties  offer  important  clues  about 
their  formation  processes  (Zinnecker  &  Yorke  2007).  However, 
our  knowledge  about  the  numbers  and  distributions  of  binary 
and  multiple  stars  is  incomplete  because  the  systems  are  gener¬ 
ally  so  distant  that  we  cannot  detect  binaries  in  the  separation 
realm  between  spectroscopically  detected  systems  and  angularly 
resolved  systems,  i.e.,  those  with  small  angular  separations  and 
periods  in  the  range  of  years  to  decades  (Mason  et  al.  1998). 
We  need  milllarcsecond  (mas)  resolution  to  start  to  fill  in  this 
observational  period  gap  (Sana  et  al.  2008). 

The  Orion  Nebula  cluster  provides  the  closest  example  (d  = 
0.41  kpc;  Menten  et  al.  2007)  of  an  environment  with  a  modest 
number  of  massive  0-stars  (Preibisch  et  al.  1999;  Weigelt 
et  al.  1999;  Close  et  al.  2012).  However,  to  explore  a  large 


Current  Address:  Department  of  Physics  and  Astronomy  University  of 
Sheffield,  Hounsfield  Road,  Sheffield  S3  7RH,  UK. 


sample  of  very  massive  stars,  the  next  closest  environment  is  the 
association  Cygnus  OB2  at  a  distance  of  <7  =  1.40  ±  0.08  kpc 
(Rygl  et  al.  2012).  The  Cyg  OB2  association  has  approximately 
2600  ±  400  members  (Rnodlseder  2000)  with  about  100  O- 
stars  within  the  central  1°  (Comeron  et  al.  2002;  Wright  et  al. 
2010),  making  it  one  of  the  largest  concentrations  of  OB 
stars  in  the  Galaxy.  It  is  home  to  some  of  the  most  massive 
(M  >  100  Mq;  Massey  &  Thompson  1991;  Rnodlseder  2000; 
Herrero  et  al.  2001)  and  intrinsically  brightest  stars  (MT  304; 
Massey  &  Thompson  1991)  known  in  our  Galaxy,  including 
two  of  the  rare  03-type  stars  (MT  417  and  MT  457;  Walborn 
et  al.  2002;  Walborn  1973,  respectively).  The  close  binary 
properties  of  the  Cyg  OB2  stars  have  been  studied  extensively  in 
the  spectroscopic  survey  of  Kiminki  &  Kobulnicky  (2012,  and 
references  therein).  A  few  wider  systems  have  been  identified 
through  high  angular  resolution  speckle  (Mason  et  al.  2009)  and 
imaging  observations  (Mafz  Apellaniz  2010). 

The  fine  guidance  sensors  (EGSs)  aboard  the  Hubble  Space 
Telescope  (HST)  provide  us  with  the  means  to  search  for 
visual  binaries  among  the  optically  faint,  massive  stars  of  Cyg 
OB2.  EGS  TRANS  mode  observations  can  resolve  systems 
with  separations  of  (y.'Ol-l",  or  14-1400  AU  (or  2800  AU  in 
the  case  of  MT  417  with  a  2"5  scan  length)  and  differential 
magnitudes  less  than  about  4  mag  (Horch  et  al.  2006;  Nelan  et  al. 
2012).  At  these  separations  we  are  sampling  systems  where  the 
components  would  evolve  independently  of  each  other  but  are 
nevertheless  important  in  understanding  massive  star  formation. 
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Due  to  the  short  length  of  the  scans  (~1")  the  sources  we  find 
are  highly  prohahle  to  be  true  companions.  We  will  discuss 
the  probability  of  chance  alignments  in  our  sample  in  a  future 
paper.  Nelan  et  al.  (2004)  used  the  FGS  instrument  to  search 
for  binaries  among  23  massive  stars  in  the  Carina  Nebula  region 
(d  ^  2.5kpc),  and  they  discovered  five  new  binaries,  including  a 
companion  of  the  very  massive  star  HD  93 129 A  at  a  separation 
of  53  mas  (133  AU).  Here  we  present  the  results  of  an  FGS 
survey  for  visual  binaries  around  58  stars  in  Cyg  OB2.  We 
compare  the  results  with  those  from  an  infrared  adaptive  optics 
survey  that  will  be  presented  in  a  future  paper.  In  Section  2 
we  describe  the  sample  selection,  observations,  and  the  FGS 
reduction  pipeline  and  process.  Section  3  describes  how  binary 
stars  are  detected,  and  Section  4  details  the  fitting  routines  used 
to  determine  the  system  parameters  (differential  magnitude, 
angular  separation,  and  position  angle).  The  detection  limits  for 
our  sample  are  also  discussed  there.  Our  results  are  presented 
in  Section  5,  and  a  discussion  of  the  multiplicity  of  our  sample 
and  future  work  is  given  in  Section  6. 

2.  SAMPLE  AND  OBSERVATIONS 

The  stars  were  selected  from  among  the  brightest,  most 
massive  stars  in  Cyg  OB2  as  cataloged  in  Schulte  (1958), 
Massey  &  Thompson  (1991),  and  Comeron  et  al.  (2002).  All  the 
stars  in  our  sample  have  published  spectral  classifications  and 
are  brighter  than  V  =  14  mag,  well  within  the  detection  limits 
of  EGS.  The  observations  were  scheduled  under  HST  proposal 
10612  (PI:  D.  Gies),  a  SNAP  program  with  70  available  targets. 
SNAP  programs  are  used  when  needed  to  fill  gaps  in  HSTs 
observing  schedule  which  cannot  be  filled  by  the  GO  class 
programs.  Typically  50%  of  the  targets  in  a  SNAP  program 
are  actually  observed;  our  program  fared  better,  with  58  stars 
observed,  listed  in  Table  1.  The  stars  are  identified  according 
to  the  numbering  scheme  used  in  the  optical  photometric 
survey  of  Massey  &  Thompson  (1991;  e.g.,  MT  138)  and  the 
infrared  spectroscopic  study  of  Comeron  et  al.  (2002;  e.g.,  A 
23).  Stars  not  included  in  those  surveys  are  identified  by  the 
number  assigned  by  Schulte  (1958;  e.g.,  SCHULTE  5  =  S  5 
=  Cyg  OB2-5)  or  by  the  Wolf-Rayet  number  from  van  der 
Hucht  (2001;  e.g.,  WR  145).  Table  1  also  lists  the  celestial 
coordinates  from  the  Two  Micron  All  Sky  Survey  (Skrutskie 
et  al.  2006),  spectral  classification  and  reference  source,  V 
magnitude,  B—V  color,  date  of  observation,  the  name  of  the 
single  star  whose  observation  was  used  as  the  calibrator  for 
the  binary  fitting,  and  the  number  of  components  detected.  The 
stars  from  Comeron  et  al.  (2002)  do  not  have  known  B—V 
colors,  but  because  they  are  bright  infrared  sources,  they  are 
assumed  to  be  redder  than  the  other  stars  in  the  survey.  The 
final  column  indicates  known  binary  systems  detected  through 
our  adaptive  optics  program  or  through  the  spectroscopic 
observations  of  Kiminki  et  al.  (2012).  “NIRI”  in  the  remarks 
column  denotes  systems  that  we  resolved  in  a  A'-band  survey 
made  with  the  NIRI  camera  and  Altair  adaptive  optics  system 
on  the  Gemini  North  Telescope  (Caballero-Nieves  2012).  “RV 
constant”  indicates  objects  observed  by  Kobulnicky  et  al.  (2012) 
that  did  not  show  radial  velocity  variability  during  the  course  of 
their  survey.  Spectroscopic  binaries  are  denoted  by  “SBl”  and 
“SB2”  for  those  detected  by  single-lined  or  double-lined  radial 
velocity  variations,  respectively  (Kiminki  et  al.  2012).  Eclipsing 
binaries  are  denoted  by  “EW/KE”  (W  UMa  type  with  ellipsoidal 
variations  and  periods  P  <  1  day),  “EA”  (detached  systems 
with  flat  maxima),  and  “EB”  (semi-detached  with  ellipsoidal 
light  curves;  Kiminki  et  al.  2012).  At  the  time  of  writing,  32 


stars  out  of  the  58  observed  are  known  multiple  systems;  9  are 
angularly  resolved  systems  listed  in  the  Washington  Double  Star 
catalog  (Mason  et  al.  2001),  two  of  which  are  also  spectroscopic 
binaries  from  the  25  systems  detected  by  Kiminki  et  al.  (2012). 

All  the  stars  were  observed  using  FGSlr  in  its  high  angular 
resolution  TRANSEER  mode  (TRANS).  The  EGS  is  a  Koesters 
prism  based,  white  light,  shearing  interferometer  that  is  sensitive 
to  the  tilt  of  the  incident  wave  front  from  a  luminous  object. 
The  EGS  optical  train  includes  a  polarizing  beam  splitter  that 
illuminates  two  mutually  orthogonal  Koesters  prisms  which 
provides  simultaneous  sensitivity  to  the  wavefront  tilt  along 
the  EGS’s  X-  and  y-axes.  In  TRANSEER  mode  the  FGS  5"  x  5" 
instantaneous  field  of  view  (IFOV)  is  scanned  across  the  object 
of  interest  along  a  45  degree  path  in  the  FGS  (x,y)  coordinate 
frame,  which  varies  the  incoming  wavefront  tilt  along  both  axes. 
If  the  source  is  a  resolved  binary,  the  light  from  the  two  stars 
will  be  mutually  incoherent  and  the  observed  interference  fringe 
will  differ  from  that  of  an  unresolved  point  source  (Nelan  et  al. 
2012).  We  observed  each  star  using  20  scans  with  a  scan  length 
of  1"  composed  of  1  mas  steps  (for  MT  417  we  employed  a 
2".  5  scan  length  to  capture  each  component  of  this  known  wide 
binary).  All  observations  were  made  using  the  F583W  filter, 
which  has  a  central  wavelength  of  5830  A  and  a  passband  of 
3400  A. 

The  resultant  interference  fringes,  traditionally  referred  to  as 
“^-curves”  due  to  their  characteristic  shape,  are  reconstructed 
from  the  photon  counts  reported  by  the  photo-multiplier  tubes 
(PMTs)  at  each  step  of  the  scan  using  the  FGS  data  reduction 
pipeline.  The  reconstructed  individual  scans  are  cross  correlated 
to  optimize  mutual  alignment  (eliminating  spacecraft  drift). 
Once  aligned,  the  scans  are  co-added  and  then  smoothed  using 
a  piecewise  spline  fit  to  obtain  the  optimal  signal-to-noise  ratio 
of  the  final  S-curve.  Any  scan  found  to  be  excessively  noisy 
(from  high  spacecraft  jitter)  was  deleted  from  the  process.  We 
found  that  it  was  helpful  to  remove  the  low  frequency,  slowly 
varying  background  of  the  scans  because  of  their  increasing 
departure  from  zero  difference  with  increasing  separation  from 
the  central  fringe  (caused  by  spatial  sensitivity  variations  of 
the  photomultiplier  tubes).  This  was  usually  accomplished  by 
subtracting  a  parabolic  fit  of  the  outer  fringe  pattern,  or  wings, 
determined  from  the  mean  of  many  calibrator  (single)  star  scans. 
However,  in  a  few  cases  the  outer  fringe  variation  was  too  large 
to  be  rectified  this  way,  and  a  spline  fit  was  made  of  the  fringe 
variation  at  representative  points  along  the  ^-curve.  The  final 
5-curves  for  all  of  the  stars  are  presented  in  alphabetical  order  in 
Figure  1  (available  in  full  in  the  electronic  version  of  the  paper). 
The  5-curves  appear  in  the  middle  panels  for  both  the  x-  and 
y-axes  while  the  upper  and  lower  panels  summarize  searches  for 
resolved  and  blended  companions,  respectively  (see  Section  3). 

3.  COMPANION  DETECTION 

Binaries  display  an  5-curve  that  appears  as  the  sum  of  two 
offset  and  rescaled  fringe  patterns  (see  Section  3.2).  A  binary 
detection  is  made  by  comparing  the  5-curve  of  a  target  to  that  of  a 
point  source  (single  star).  We  describe  below  the  three  detection 
methods  we  used  to  determine  binarity.  The  comparison  star,  or 
calibrator,  is  taken  from  unresolved  sources  within  our  sample 
that  meet  three  criteria:  (1)  observed  using  the  same  filter,  (2) 
close  'mB—V  color,  and  (3)  close  in  observation  date.  The  first 
criterion  is  met  by  all  the  stars  in  our  sample.  Criterion  2  is 
necessary  because  the  appearance  of  a  point  source  5-curve  is 
color  dependent  due  to  refractive  optics  within  the  EGS  and 
the  wavelength  dependent  response  of  the  PMTs.  Eor  example. 
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Table  1 

Cygnus  OB2  Target  List 


Caballero-Nieves  et  al. 


Star 

Name 

(1) 

R.A. 

(J2000) 

(2) 

Deck 

(J2000) 

(3) 

Spectral 

Classification 

(4) 

Class. 

Ref. 

(5) 

V 

(mag) 

(6) 

B-V 

(mag) 

(7) 

Obs.  Date 
(BY) 

(8) 

Calibrator 

Name 

(9) 

No. 

Comp. 

(10) 

Remarks 

(11) 

A  23 

20:30:39.71 

+41:08:49.0 

B0.7  Ib 

1 

11.25 

2007.2729 

1 

A  27 

20:34:44.72 

+40:51:46.6 

BO  la 

2 

11.26 

2007.4675 

1 

A41 

20:31:08.38 

+42:02:42.3 

09.7  II 

1 

11.70 

2006.0666 

2 

NlRl 

A  46 

20:31:00.20 

+40:49:49.7 

07  V((f)) 

1 

11.40 

2008.4988 

1 

MT5 

20:30:39.82 

+41:36:50.7 

06  V 

3 

12.93 

1.64 

2006.5055 

WR  145 

2 

NlRl 

MT59 

20:31:10.55 

+41:31:53.5 

08  V 

4 

11.18 

1.47 

2006.4855 

MT601 

2 

NIRl,  SBl 

MT70 

20:31:18.33 

+41:21:21.7 

09  II 

4 

12.99 

2.10 

2006.4928 

1 

SBl 

MT83 

20:31:22.04 

+41:31:28.4 

B1  I 

3 

10.64 

1.18 

2006.0001 

1 

SBl 

MT  138 

20:31:45.40 

+41:18:26.8 

081 

3 

12.26 

1.99 

2006.4854 

MT390 

2 

NIRl,  SBl 

MT  145 

20:31:49.66 

+41:28:26.5 

09  III 

4 

11.52 

1.11 

2005.9990 

1 

SBl 

MT213 

20:32:13.13 

+41:27:24.6 

BOV 

3 

11.95 

1.13 

2007.2387 

1 

RV  constant 

MT217 

20:32:13.83 

+41:27:12.0 

07  Illf 

3 

10.23 

1.19 

2006.3752 

1 

RV  constant 

MT227 

20:32:16.56 

+41:25:35.7 

09  V 

3 

11.47 

1.24 

2006.4883 

1 

MT250 

20:32:26.08 

+41:29:39.4 

B2III 

3 

12.88 

1.06 

2007.2942 

1 

RV  constant 

MT258 

20:32:27.66 

+41:26:22.1 

08  V 

4 

11.10 

1.20 
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Figure  1.  Final  5-curves  for  the  x  and  y  orthogonal  scans  for  A  23.  The  top  panels  show  the  mean  residual  cross-correlation  functions  (black  line:  after  removal  of 
the  primaiy  component)  plus  the  CCF  standard  deviations  (gray  line).  The  positions  of  components  are  indicated  by  vertical  dashed  lines.  The  middle  panels  show 
the  observed  (solid  line)  and  modeled  (dashed  line)  5-curves.  The  lower  panels  display  the  observed  minus  model  residuals  (black  line)  and  standard  deviation  (gray 
region)  plus  a  ht  of  the  second  derivative  (gray  line)  if  a  blended  companion  is  present. 

(The  complete  figure  set  (58  images)  is  available  in  the  online  journal.) 


the  shape  of  an  5-curve  is  slightly  broadened  for  redder  targets 
(Horch  et  al.  2006).  We  chose  calibrators  whose  colors  are  within 
±0.2  mag  from  the  target  B  —V,  except  for  a  few  very  red  stars 
(Section  4).  The  appearance  of  the  5-curves  can  change  with 
time  due  to  the  settling  of  the  instruments  and  adjustments  from 
servicing  missions.  Our  observing  program  was  not  affected 
by  any  changes  made  due  to  a  servicing  mission,  but  our  large 
range  of  observation  dates  (2005  December  to  2008  June)  spans 
a  range  where  the  long  term  changes  in  FGSlr,  while  small, 
are  non-negligible.  Thus,  we  chose  as  the  best  calibrator  the  one 
observed  closest  in  time  to  the  target.  The  binary  systems  we  are 
able  to  resolve  have  components  with  nearly  equal  brightness 
and  hence,  their  colors  should  be  comparable  in  most  cases. 
Consequently,  we  adopted  the  same  calibrator  to  model  both 
components. 

3.1.  Visual  Inspection 

The  5-curve  or  transfer  function  of  a  binary  consists  of  the 
normalized  superposition  of  the  individual  transfer  functions 
of  two  point  sources.  For  widely  separated  systems,  the  scan 


will  show  two  shifted  transfer  functions,  whose  relative  am¬ 
plitudes  are  directly  related  to  the  magnitude  difference  (see 
Equation  (3)).  For  closer  systems,  the  5-curve  can  look  obvi¬ 
ously  different  from  that  of  a  single  star  (such  as  the  case  of 
MT  429  shown  in  Figure  1.25),  because  the  composite  fringe  is 
formed  from  two  overlapping  fringes.  Thus,  a  direct  comparison 
of  the  smoothed,  co-aligned,  target  5-curve  with  that  of  a  single 
star  provides  the  first  way  to  identify  binaries. 

Another  indicator  of  a  resolved  binary  is  the  relative  fringe 
amplitude,  calculated  as  the  value  of  sppr,  or  the  5-curve  peak- 
to-peak  ratio.  The  ratio  of  the  extrema  points  of  the  fringe  to 
that  of  a  single  star  calibrator  fringe  can  be  expressed  as 

•^obs.max  ‘^obs.min 

sppr  =  — (1) 
.max  —  5'f583W  ,min 

where,  for  bright  stars  (V  <  12),  5F583w,max  =  0.61  and 
>S'F583w,min  =  -0.55  along  the  x-axis  and  5F583w,max  =  0.30 
and  5F583w,min  =  —0.65  along  the  y-axis  for  typical  single 
star  observations  made  with  the  F583W  filter.  This  parameter 
is  listed  in  the  middle  panels  of  Figure  1  using  denominator 
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values  from  the  mean  of  the  selected  calibrator  5-curves. 
The  cases  where  sppr  ^  0.9  indicate  that  the  target  is  a 
resolved  binary  in  which  destructive  interference  causes  a 
decline  in  the  overall  fringe  amplitude.  However,  if  0.9  ^ 
sppr  ^  1.0,  as  would  be  the  case  for  a  binary  with  a 
small  angular  separation  or  large  magnitude  difference  of 
the  components,  a  more  thorough  analysis  is  needed.  Note 
that  the  x-  and  y-axes  have  different  values;  this  is  due  to 
the  effect  of  the  HST  spherical  aberration  and  the  alignment 
of  the  interferometer  optical  elements  relative  to  the  HST 
optical  axis.  Also  note  that  the  5f853w  maximum  and  minimum 
values  are  adjusted  for  the  magnitude  of  the  target  since  the 
5-curve  is  also  affected  by  the  PMT  dark  current,  which  becomes 
an  increasingly  larger  percentage  of  the  total  counts  as  fainter 
targets  are  observed.  This  effect  is  negligible  for  V  <  12,  but  is 
significant  for  fainter  objects.  Note,  the  initial  sppr  value  used 
as  a  quick-look  tool  for  a  companion  was  computed  from 
the  scans  that  have  not  been  de-jittered,  cross  correlated,  or 
smoothed  as  they  are  for  the  more  thorough  analysis  described 
in  the  next  sections.  The  sppr  values  quoted  in  Figure  1  are 
calculated  after  this  analysis  has  been  performed,  but  do  not 
differ  significantly  from  the  initial  sppr  value  used. 

We  produced  an  initial  list  of  stars  that  appeared  single 
according  to  the  visual  inspection  and  the  sppr  filter  as  possible 
calibrators  for  the  analysis  described  in  the  following  sections. 
From  there,  we  narrowed  the  list  through  an  iterative  process  of 
selecting  which  stars  met  the  following  criteria  for  point  sources, 
and  we  used  those  as  the  calibrators  for  modeling  purposes. 

3.2.  Marginally  Resolved  Systems  and  Detection  Limits 

Close  systems  (projected  separation  <25  mas),  where  the  two 
5-curves  are  not  well  separated,  can  slightly  reduce  the  fringe 
amplitude  and  widen  the  fringe  shape.  If  the  fringe  pattern  of 
a  single  calibrator  star  is  5(x),  then  the  observed  fringe  pattern 
for  more  than  one  star  will  be 

n 

S(x)obs  Six  -  Xi)  (2) 

where  each  of  n  stars  has  a  flux  fraction  /)  =  Fj  and 

a  relative  projected  offset  position  x,-.  For  a  binary  star  with  a 
companion  flux  ratio  r  —  F^jFi  and  a  projected  separation  Ax, 
the  observed  pattern  simplifies  to 

5(x)obs  =  - 5(x)  H-  - - 5(x  -  Ax).  (3) 

\+r  1+r 

An  analytical  representation  of  the  difference  between  the 

binary  and  calibrator  5-curves  can  be  estimated  by  making  a 
second-order  expansion  for  small  offset  e, 

5(xh-6)  =  5(x)H-e5'(x)H-^e^5"(x)  (4) 

where  S'  and  5"  are  the  hrst  and  second  derivatives  of  the 
5-curve.  In  the  frame  of  reference  where  5(0)  =  0  for  the  binary, 
the  primary  and  secondary  5-curves  will  be  respectively  shifted 
by  amounts  ei  =  — (r/(l  H- r))Ax  and  €2  —  H-(l/(l  H- r))Ax, 
where  Ax  is  the  projected  separation  of  secondary  from  primary. 
Then  the  difference  between  the  marginally  resolved  binary  and 
calibrator  5-curves  is 

1  r 

5(x)bin  -  5(x)cai  =  T - Six  +  ei)  H-  - - 5(x  h-  €2)  -  Six) 

1  +  r  \+r 

=  5  (5) 

2(1  +r)2 


X  (arcsec) 


Figure  2.  Example  plots  of  the  difference  between  a  model  binary  and  calibrator 
^-curves.  The  dashed  line  shows  the  difference  for  a  binaiy  with  r  =  I  and 
Ax  =  (y/015,  and  the  dot-dashed  line  shows  the  difference  for  a  binaiy  with 
r  =  0.5  and  Ax  =  0'.'0159.  Both  resemble  the  second  derivative  of  the  calibrator 
^-curve  (solid  line)  for  the  same  amplitude  a,  but  differ  from  that  caused  by  the 
dilution  of  an  off-scan  companion  (dotted  line). 

This  second-order  expression  has  several  important  features. 
First,  the  observed  difference  in  the  core  of  the  5-curve  will 
appear  to  have  the  same  functional  shape  as  the  second  derivative 
of  the  5-curve,  so  we  can  directly  search  for  companions  with 
overlapping  fringes  by  looking  for  a  difference  that  has  a 
second  derivative  shape.  Second,  the  amplitude  of  the  difference 
depends  on  a  product  involving  both  the  separation  Ax  and  the 
flux  ratio  r,  so  in  the  absence  of  other  information,  neither 
parameter  can  be  determined  uniquely.  Third,  the  amplitude 
of  the  difference  depends  on  the  separation  squared,  so  no 
information  can  be  reliably  extracted  on  the  direction  of  the 
companion  from  the  primary. 

Figure  2  shows  examples  of  such  5-curve  differences  for 
model  binaries.  The  dashed  line  shows  the  difference  for  a  model 
of  equally  bright  stars  (r  =  1)  with  a  separation  of  (f.'OlS  made 
from  a  mean  5-curve  from  a  collection  of  calibrator  x-axis  scans. 
According  to  the  analytical  expression  above,  the  coefficient 
leading  the  second  derivative  is  0.015^/8  =  2.8  x  10“^  arcsec^, 
and  the  solid  line  shows  the  product  of  this  coefficient  and 
a  numerical  solution  of  the  second  derivative  of  the  cali¬ 
brator  5-curve  (smoothed  by  convolution  with  a  Gaussian  of 
FWHM  =  O'.' 005).  The  good  match  between  the  detailed  model 
and  analytical  solution  verifies  the  second  derivative  charac¬ 
ter  of  the  difference  curve.  The  same  coefficient  is  found  for 
r  =  0.5  and  Ax  =  0'.'0159,  and  the  dot-dashed  line  shows  the 
difference  of  the  binary  and  calibrator  curves  for  these  binary 
parameters.  Again,  the  agreement  between  this  model  and  the 
analytical  curve  shows  that  two  models  with  the  same  product 
a  =  (l/2)(r/(l  +  r)^)(Ax)^  have  very  similar  5-curves. 

Therefore,  in  the  case  of  marginally  resolved  systems  the 
difference  between  the  5-curves  of  a  suspected  binary  and  a 
single  star  should  look  like  that  of  the  second  derivative  of  the 
point-source  transfer  function  scaled  by  the  coefficient  product 
term  a.  Unless  the  flux  ratio  is  determined  independently,  there 
is  not  a  unique  solution  for  r  and  Ax.  The  method  was  applied 
by  considering  the  difference  between  the  target  and  calibrator 
5-curves  over  the  range  within  ±100  mas  of  the  center  of  the 
fringe.  The  coefficient  was  then  estimated  by  a  least-squares 
fit  of  Equation  (5)  over  the  restricted  range.  The  coefficient 
a  was  determined  in  practice  with  an  ensemble  of  like-color 
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Figure  3.  Example  of  the  cross-correlation  function.  The  top  panel  shows  an 
example  of  a  model  S-curve  of  a  binary  star  with  p  =  T-O'.'OV  and  a  flux  ratio 
r  =  F2IF1  =  0.1  (solid  line)  and  that  of  a  single  star  calibrator  (dotted  line).  The 
lower  panel  shows  the  cross-correlation  functions  of  the  target  with  the  calibrator 
(solid  line)  and  of  the  calibrator  with  itself  (dotted  line).  The  difference  between 
these  two  is  shown  as  a  dashed  line  in  an  expanded  scale,  offset  for  clarity. 

calibrators,  and  the  criterion  for  detection  was  set  by  a  mean 
coefficient  a  with  a  positive  value  greater  than  4cr ,  where  a  is 
the  standard  deviation  of  the  coefficient  derived  using  different 
calibrators. 

3.3.  Wide  Binary  Detection 

The  next  situation  to  consider  is  the  case  where  the  absolute 
projected  separation  of  the  binary  companion  is  comparable  to 
or  greater  than  the  width  of  the  fringe  (ss50  mas  or  70  AU  for 
stars  at  the  distance  of  Cyg  OB2)  and  the  secondary  may  be 
faint.  The  best  approach  in  these  cases  is  to  calculate  the  cross¬ 
correlation  function  (CCF)  of  a  target  ^-curve  with  that  of  a 
calibrator  star.  This  method  has  the  advantage  of  using  more  of 
the  ^-curve  than  just  the  extrema  points,  and  it  potentially  helps 
unravel  those  cases  where  the  fringes  overlap.  The  cost  of  this 
approach  is  a  slight  decrease  in  the  working  angular  resolution 
limit  (but  see  above  for  a  discussion  of  binaries  with  blended 
^-curves).  The  top  panel  of  Figure  3  shows  an  example  of  a 
model  S'-curve  of  a  binary  star  with  a  projected  separation  of 


-t-fy.'O?  and  a  flux  ratio  r  —  F^jFi  —  0.1  (constructed  using 
calibrator  x-axis  scans).  The  dotted  line  shows  the  calibrator 
5-curve  while  the  solid  line  shows  that  for  the  target  binary.  The 
fringe  patterns  overlap  significantly  at  this  separation,  and  the 
main  differences  are  a  dilution  of  the  main  fringe  pattern  and  a 
change  in  outer  fringe  structure  near  X  =  H-O'.'OS.Thelowerpanel 
shows  the  CCFs  of  the  target  with  the  calibrator  (solid  line)  and 
of  the  calibrator  with  itself  (dotted  line).  Unlike  the  5-curves,  the 
CCFs  show  one  main  peak  for  each  stellar  component.  In  order 
to  isolate  the  companion,  the  calibrator  CCF  is  shifted  to  the 
peak  of  the  target  CCF  and  scaled  to  the  peak  of  the  target  CCF. 
The  shifted  and  scaled  calibrator  CCF  is  then  subtracted  from 
the  target  CCF  to  produce  the  residual  CCF  shown  as  a  dashed 
line  in  an  expanded  scale  in  the  lower  panel  (and  offset  by  —0.8 
for  clarity).  Now  the  peak  from  the  companion  is  clearly  visible 
at  the  offset  position  of  Ax  =  h-OTO?. 

We  need  a  working  criterion  to  establish  whether  or  not 
a  peak  in  the  residual  CCF  makes  a  significant  detection  of 
a  companion.  Because  the  dominant  source  of  uncertainty 
in  the  shape  of  the  5-curves  is  the  inherent  scatter  between 
observations  of  the  calibrators,  the  criterion  was  set  by  running 
the  CCF  procedure  for  any  given  target  with  an  ensemble  of 
calibrator  5-curves  for  stars  of  similar  color  (usually  a  set  of  ten 
calibrators).  Then  the  detection  criterion  was  set  by  requiring 
the  peak  in  the  mean  of  the  residual  CCFs  to  exceed  4ct(x), 
where  (t(x)  is  the  standard  deviation  of  the  residual  CCFs  at  the 
peak  position  x. 

Figure  1  shows  transfer  function  plots  for  all  58  stars  in  our 
sample.  The  hgures  show  the  x-  (left)  and  y-axis  (right)  rectified 
5-curves  in  the  central  panel.  Also  shown  as  a  dashed  line  is  a 
preliminary  model  ht  based  upon  the  components  derived  from 
the  CCF  analysis  and  the  mean  5-curve  of  the  calibrator  set 
selected.  The  top  panel  plots  the  mean  of  the  residual  CCFs  of 
the  target  with  the  calibrator  (solid  black  lines)  after  the  peak  of 
the  primary  has  been  subtracted  off.  The  vertical  dashed  lines 
indicate  the  position  of  each  component  resolved.  The  solid, 
gray  lines  show  the  standard  deviation  of  the  residual  CCFs,  and 
a  peak  is  considered  significant  only  if  the  mean  residual  CCF 
(black)  exceeds  the  standard  deviation  (gray)  by  four  times.  The 
bottom  panels  show  the  difference  curves  between  the  star  and 
model  calibrator  5-curves  (solid  line).  The  shaded,  gray  region  is 
the  uncertainty  envelope  determined  from  the  standard  deviation 
at  each  point  along  the  5-curves  for  the  calibrators.  In  the  case 
where  the  second  derivative  test  resolved  a  blended  pair  (e.g., 
MT  304,  Figure  1.18),  the  second  derivative  of  the  calibrator 
5-curve  is  overplotted  and  scaled  by  the  a  coefficient  (solid 
gray  line). 

The  detection  threshold  for  a  binary  of  a  given  separation  is 
determined  by  comparing  multiple  examples  of  a  model  binary 
based  upon  the  calibrator  observations  to  the  calibrator  5-curves. 
This  was  done  using  a  set  of  21  calibrator  observations  of  similar 
color  stars  for  both  the  model  and  calibrator  curves  (for  a  total 
of  420  test  cases),  and  the  faintest  flux  ratio  was  set  by  models 
that  met  the  4cr  detection  criterion.  The  positive  (solid  line)  and 
negative  (dashed  line)  branches  are  folded  onto  one  separation 
axis  in  Figure  4  in  the  left  and  right  panels  for  the  x-  and  y-axis 
scans,  respectively.  Here  the  limiting  flux  ratio,  r,  is  shown  as 
a  magnitude  difference  Am  =  —2.5  log  r  for  trial  separations, 
and  any  binary  brighter  than  the  limit  (i.e.,  below  the  line  plots) 
would  exceed  the  4a  detection  criterion.  The  smallest  separation 
detected  in  favorable,  equal  flux  cases  is  slightly  better  for 
positive  separations,  and  the  faintest  detectable  companions 
have  Am  =  4.5  mag  at  large  separations.  The  oscillation  in 
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jr  PROJECTED  SEPARATION  (arcsec)  y  PROJECTED  SEPARATION  (arcsec) 

Figure  4.  Binary  detection  limits  as  a  function  of  projected  separation  and  magnitude  difference.  The  solid  (dashed)  lines  indicate  the  limits  for  positive  (negative) 
offsets  in  x  (left  panel)  and  y  (right  panel),  respectively.  Binaries  with  separations  and  magnitude  differences  below  these  limits  should  be  detected.  The  dotted  lines 
indicate  detection  limits  for  marginally  resolved  binaries  using  the  second  derivative  method. 


these  curves  seen  near  x  =  y  =  (y!06  is  due  to  the  changing 
and  relatively  larger  uncertainties  in  the  calibrator  S'-curves  at 
such  distances  from  their  zero-crossing.  The  results  in  these 
figures  compare  well  with  the  advertised  limits  in  Figure  3.3  of 
the  Fine  Guidance  Sensor  Instrument  Handbook  (Nelan  et  al. 
2012)  except  for  the  case  of  marginally  resolved  binaries  that 
we  discussed  above. 

Note  that  the  CCF  method  loses  its  effectiveness  for  very 
close  companions.  The  peaks  in  the  residual  CCFs  in  the  binary 
models  where  p  <  25  mas  are  rarely  found  at  separations 
of  less  than  O'.' 04  even  in  favorable  cases.  This  is  due  to  the 
fact  that  blending  of  the  fringe  patterns  becomes  so  severe 
that  the  calibrator  CCF  is  positioned  at  the  maximum  that 
occurs  between  the  actual  positions  of  the  components,  and 
consequently  the  residual  CCF  shows  two  peaks:  one  for  the 
companion  and  a  mirror  one  for  the  primary.  At  the  smallest 
separations  where  the  method  can  be  applied  (WOlb  in  x  and 
(y.'Ol?  in  y),  the  two  peaks  in  the  residual  CCF  approach  equal 
intensity,  and  the  method  can  no  longer  distinguish  the  direction 
of  primary  to  secondary.  Nevertheless,  the  appearance  of  a 
double-peaked  feature  in  the  residual  CCF  offers  evidence 
for  the  presence  of  a  marginally  resolved  companion  that 
supplements  the  second  derivative  test. 

The  detection  limits  using  the  second  derivative  method 
were  estimated  by  running  the  scheme  with  multiple  binary 
models  from  a  large  sample  of  calibrators.  The  working  criteria 
from  these  models  for  4a  detection  are  a  >  1.2  x  10  ^  and 
1.6  X  10“^  arcsec^  for  the  x-  and  y-axes,  respectively.  These 
limits  are  shown  in  Figure  4  as  dotted  lines  that  trace  the  upper 
envelope  for  detection  by  the  second  derivative  approach.  In 
the  best  cases  (r  =  l),theestimatessuggestthatbinariesasclose 
as  O'.'OlO  can  be  detected  with  FGS  TRANS  mode  scans.  The 
second  derivative  and  CCF  methods  are  probably  both  sensitive 
to  binary  detection  in  the  O'.' 020-0'.' 025  (28-100  AU)  range. 

3.4.  Off-scan  Components 

There  is  also  the  case  of  very  widely  separated  binaries,  where 
the  light  of  the  companion  falls  within  the  instrumental  FOV 
(5"  X  5"),  but  the  projected  separation  is  greater  than  the  length 
of  the  scan.  This  is  the  case  for  our  observation  of  MT  531  where 
the  system  is  resolved  along  the  y-axis,  but  the  secondary  is  off 


the  scan  along  the  x-axis  (see  Figure  1.39).  The  CCF  method 
will  fail  to  detect  the  companion  because  its  peak  lies  beyond 
the  recorded  scan.  However,  such  binaries  will  still  cause  the 
5-curve  of  the  primary  to  appear  with  an  amplitude  reduced 
by  a  factor  1/(1  H-  r),  and  hence  the  ratio  of  target  to  calibrator 
5-curve  amplitude  (sppr)  provides  an  additional  criterion  to 
check  for  very  wide  binaries.  In  practice,  the  scatter  between 
calibrator  5-curves  of  similar  color  indicates  that  a  4a  detection 
may  be  claimed  if  the  target  5-curve  amplitude  is  less  than 
92%  of  the  mean  calibrator  5-curve  amplitude.  Recall  that 
overlapping  fringes  due  to  a  binary  may  also  cause  the  amplitude 
to  decline,  but  in  this  case,  the  fringe  will  also  be  widened. 
Consequently,  one  can  differentiate  between  the  off  scan  and 
blended  cases  by  the  appearance  of  the  difference  curve.  The 
dotted  line  in  Figure  2  shows  that  the  difference  curve  for  a 
binary  with  an  off-scan  companion  appears  like  a  negatively 
scaled  version  of  the  5-curve  itself,  which  looks  very  different 
from  the  second  derivative  curve.  Thus,  it  is  important  to  inspect 
the  shapes  of  the  difference  curves  in  order  to  decide  if  a  positive 
detection  indicates  the  presence  of  an  off-scan  companion  or  a 
very  close  companion. 

We  were  able  to  make  a  second  check  for  wide  companions 
through  inspection  of  results  from  our  adaptive  optics  (AO) 
survey  using  the  Near  InfraRed  Imager  and  Spectrograph  (NIRI) 
at  the  Gemini  North  Observatory  (Caballero-Nieves  2012).  The 
NIRI  observations  have  a  larger  dynamic  range  (i.e.,  can  detect 
fainter  companions)  and  have  a  larger  field  of  view,  so  any 
system  observed  with  FGS  with  p  >  50  mas  will  also  be 
detected  in  the  NIRI  image.  The  infrared  images  allow  us  to 
identify  those  companions  with  projected  separations  longer 
than  the  scan  length,  but  still  within  the  FOV  of  the  FGS.  With 
the  help  of  the  NIRI  observations,  we  were  able  to  conclude  that 
distant  companions  influence  the  FGS  results  for  MT  59,  MT 
138,  and  MT531. 

4.  MODEL  FITTING 

If  a  binary  was  detected,  then  we  fit  the  target  5-curve  with 
a  model  binary  formed  from  a  calibrator  5-curve  in  order 
to  determine  the  relative  brightness  and  projected  separation 
along  an  axis.  We  used  two  routines  to  calculate  the  best 
fit  model.  The  first  is  BINARY_FIT  which  is  part  of  the 
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STScl  reduction  package  based  on  algorithms  developed  by 
the  Space  Telescope  Astrometry  Team  at  Lowell  Observatory 
(Franz  et  al.  1991).  In  addition,  we  developed  the  interactive 
data  language  (IDL)  routine  TRIPFIT  for  the  special  cases 
where  BINARY_FIT  was  not  able  to  converge  to  a  satisfactory 
solution.  Both  BINARY_FIT  and  TRIPFIT  use  the  x  and  y 
projected  separations  and  the  telescope  roll  angle  to  determine 
the  binary  separation  p  and  position  angle  9  (measured  east 
from  north). 

Each  binary  was  compared  to  models  based  on  four  calibra¬ 
tors.  The  calibrators  were  selected  to  have  a  B—V  color  within 
±0.2  mag  of  the  target’s  color.  The  only  exception  was  MT  304, 
with  a  B  —  V  —  3.35,  which  is  the  reddest  star  in  our  sample 
by  more  than  1  mag.  In  this  case  we  selected  calibrators  from 
the  reddest  single  star  in  our  sample,  MT  448,  and  the  stars 
from  Comerdn  et  al.  (2002;  A  23,  A  27,  A  41,  A  46),  which 
are  red  objects  according  to  their  brightness  in  the  infrared. 
We  adopted  the  fit  made  using  the  calibrator  closest  in  color 
and  observation  date  to  determine  the  system  parameters.  The 
spread  in  results  from  fits  made  using  the  other  calibrators  was 
used  to  determine  the  fitting  parameter  uncertainties.  The  same 
calibrator  was  used  to  model  both  the  primary  and  secondary 
j'-curves,  i.e.,  we  assumed  that  any  color  difference  between  the 
components  is  negligible.  Our  calibrator  selections  are  indicated 
in  column  9  of  Table  1 . 

The  program  BINARY_FIT  uses  a  least  squares  approach  to 
determine  the  projected  separation  and  magnitude  difference  of 
the  system  (Nelan  &  Makidon  2011).  BINARY_FIT  fits  scans 
from  one  axis  at  a  time,  starting  from  initial  estimates  for  the  pro¬ 
jected  separation  along  the  axis  and  the  differential  magnitude. 
If  the  results  for  differential  magnitude  from  the  x-  and  y-axis 
solutions  agreed  within  0.2  mag,  then  the  individual  separation 
results  were  adopted  as  fit,  but  we  report  the  magnitude  dif¬ 
ference  for  the  axis  solution  with  the  larger  separation  (usually 
more  reliable).  On  the  other  hand,  if  the  differential  magnitudes 
differed  by  more  than  0.2  mag,  then  we  adopted  the  value  from 
the  axis  solution  with  the  greater  separation,  and  then  re-ht 
the  scan  for  the  other  axis  by  setting  the  adopted  magnitude 
difference. 

We  encountered  several  cases  where  BINARY_FIT  could  not 
be  used.  (1)  The  fitting  code  is  limited  by  the  scan  length  of 
the  calibrator.  If  the  separation  of  the  binary  is  larger  than 
the  scan  length  of  the  calibrator,  the  program  is  not  able  to 
recreate  a  binary  wide  enough  to  model  properly  the  target.  (2) 
BINARY_FIT  only  considers  solutions  where  both  components 
are  within  the  scan  length.  For  example,  if  the  companion  is 
recorded  on  one  axis  but  lies  off  the  scan  of  the  other,  then  the 
program  will  converge  to  a  solution  of  a  very  close  system  for 
the  axis  where  the  star  is  absent.  (3)  BINARY_FIT  cannot  be 
applied  to  systems  with  more  than  two  components,  such  as  the 
triple  MT417. 

We  coded  the  IDL  program  TRIPFIT  to  hnd  the  best  ht 
model  for  triple  systems  using  a  Levenberg-Marquardt  least- 
squares  method.  The  program  is  capable  of  making  hts  of 
triple  systems,  binaries,  and  off-scan  components.  The  user 
selects  initial  estimates  of  the  positions  of  the  components 
and  the  differential  magnitudes.  The  program  then  fits  the 
j'-curves  one  axis  at  a  time  and  returns  the  best  fit  for  the 
differential  magnitudes  and  separations.  Both  BINARY_FIT 
and  TRIPFIT  return  similar  solutions  when  modeling  a  system. 
However,  for  most  cases  we  adopted  the  results  from  the 
BINARY_FIT  models  because  it  uses  the  information  from 
both  axes  simultaneously  in  determining  the  best  fit.  The  cases 


where  we  could  not  converge  on  a  satisfactory  solution  with 
BINARY_FIT,  we  report  the  results  using  TRIPFIT  and  those 
systems  are  noted  in  Section  5. 

There  are  three  sources  of  uncertainty  in  the  analysis  of  FGS 
data.  ( 1 )  Internal  errors  arise  from  photometric  shot  noise,  which 
is  important  for  stars  with  V  >  14,  and  spacecraft  jitter  that 
cannot  be  removed  using  the  guide  star  centroids  reported  by 
the  guiding  FGSs.  (2)  The  S-curves  slowly  evolve  over  time, 
which  is  due  to  small  changes  in  the  alignment  of  the  Koesters 
prisms  relative  to  the  optical  axis  of  HST,  made  sensitive  due 
to  the  spherical  aberration  of  the  HST  primary  mirror.  By 
choosing  calibrator  stars  that  are  observed  close  in  time  (on 
the  order  of  1  yr)  to  the  science  observations  this  evolutionary 
effect  is  mitigated.  (3)  Systematic  differences  between  the 
calibrators  exist.  To  investigate  how  the  photometric  noise  and 
jitter  influence  the  derived  parameters  of  a  binary,  we  selected 
four  binary  systems  (MT  5,  MT  429,  MT  605,  and  MT  632) 
and  binned  their  scans  into  three  or  four  independent  subsets. 
Each  subset  includes  approximately  five  scans  that  were  shifted, 
co-added,  and  smoothed  using  the  same  approach  applied  to  the 
complete  set  of  scans.  The  resultant  subset  of  co-aligned  scans 
was  then  fit  with  the  best  calibrator  using  BINARY_EIT.  The 
small  spread  in  the  values  of  separation  and  magnitude  shows 
that  the  internal  error  is  not  a  significant  source  of  uncertainty  for 
the  binary  parameters  (except  in  the  case  of  MT  5;  see  Section  5). 
Consequently,  the  binary  parameter  uncertainties  are  dominated 
by  the  differences  between  the  calibrator  scans,  and  in  Table  2 
we  report  uncertainties  based  on  the  standard  deviation  of  the 
parameter  fits  made  with  the  different  calibrator  5-curves. 

5.  MULTIPLICITY  RESULTS 

Table  2  lists  the  model  fitting  results  for  separation  and 
differential  magnitude  for  the  13  resolved  systems  in  our 
Cyg  OB2  sample.  The  cited  errors  are  from  a  comparison  of 
results  from  fits  made  with  different  calibrators  as  described 
in  Section  4.  There  were  two  additional  cases  where  only  the 
second  derivative  analysis  indicated  a  possible  close  binary 
system  that  was  partially  resolved  on  only  one  axis.  The  derived 
fl-coefficient  was  a  —  (15.2  ±  3.4)  x  10“®  arcsec^  along  the 
V-axis  for  MT  227  and  a  —  (25.5  ±  3.6)  x  10“®  arcsec^  along 
the  x-axis  for  MT  317.  These  two  marginally  resolved  systems 
are  objects  of  interest  for  follow-up  analysis.  We  describe  below 
the  individual  cases  for  the  fully  resolved  systems. 

A  41.  The  CCF  analysis  of  A  41  (see  Figure  1.3)  re¬ 
veals  the  presence  of  a  faint  companion  in  the  y-axis  scan. 
This  companion  was  also  found  in  the  A'-band  NIRI  results 
at  a  separation  of  p  =  0'.'35,  essentially  the  same  as  the 
y-axis  projected  separation.  This  implies  that  the  companion’s 
projected  separation  along  the  x-axis  is  within  the  fringe  of  the 
primary.  This,  along  with  its  relative  faintness  (5  F  Rs  3.3),  make 
detection  along  the  x-axis  challenging.  Nevertheless,  we  fit  the 
5-curves  of  both  axes  with  BINARY_FIT  to  obtain  a  separation 
and  position  angle  consistent  with  the  NIRI  results. 

MT  5.  (see  Figure  1.5)  After  splitting  the  5-curve  data  into 
subsets,  the  internal  error  was  found  to  be  a  non-negligible 
source  of  uncertainty  in  this  case.  The  uncertainties  listed  for 
MT  5  in  Table  2  reflect  both  the  spread  among  the  calibrators 
and  between  the  subsets  of  MT  5  (added  in  quadrature).  Note 
that  the  second  derivative  test  for  the  x-axis  (a  =  7.9ct)  indicates 
that  the  primary  itself  may  have  a  close  companion  making  this 
a  triple  system,  but  the  test  is  consistent  with  a  single  primary 
for  the  y-axis  scan. 
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Table  2 

Multiplicity  Parameters  for  Resolved  Systems 


Star 

Name 

A.T 

(mas) 

Ay 

(mas) 

P 

(mas) 

e 

n 

AniF583W 

(mag) 

A  41 

54.4  ±  15.2 

352.6  ±  1.4 

356.8  ±  3.1 

276.1  ±  6.5 

3.27  ±  0.86 

MT5 

-129.6  ±  1.4 

294.0  ±  1.1 

321.3  ±0.8 

91.7  ±  0.3 

2.79  ±  0.09 

MT59 

-751.3  ±0.8 

>751.3  ±0.8 

2.58  ±  0.04 

MT  138 

66.0  ±3.1 

>66.0  ±3.1 

2.79  ±  0.24 

MT304 

61.8  ±  3.1 

(-14.9  or  9.7)  ±4.6 

63.6  ±  3.5 

(305.9  or  283.5)  ±  3.3 

2.31  ±0.21 

MT417  A,Ba 

-231.3  ±  1.3 

1505.0  ±  3.9 

1522.7  ±  3.7 

146.9  ±  0.5 

0.45  ±  0.01 

MT  417  A,Bb 

-214.1  ±  5.1 

1700.0  ±  2.5 

1715.7  ±2.9 

147.9  ±  0.4 

2.66  ±  0.32 

MT429 

-36.2  ±  0.7 

-95.8  ±  0.3 

101.9  ±0.3 

23.5  ±  0.4 

1.09  ±0.02 

MT516 

374.5  ±  0.04 

-618.0  ±0.1 

722.6  ±  0.09 

325.7  ±  0.1 

0.28  ±  0.02 

MT531 

470.8  ±  0.08 

>  470.8  ±  0.08 

0.50  ±  0.04 

MT605 

99.0  ±0.15 

61.2  ±0.5 

116.4  ±0.4 

255.0  ±  0.2 

0.69  ±  0.02 

MT632 

165.0  ±  0.3 

-144.1  ±0.3 

219.1  ±0.2 

247.2  ±  0.09 

2.00  ±0.16 

MT696 

(±9.7)±  1.5 

20.6  ±  2.7 

22.8  ±  3.0 

(175.1  or  225.6)  ±  1.9 

0.94  ±  0.40 

SCHULTE  5 

-801.1  ±2.2 

439.6  ±2.1 

913.9  ±2.3 

56.4  ±0.11 

2.93  ±  0.29 

MT  59.  TRIPFIT  was  used  to  model  MT  59  (see  Figure  1 .6) 
because  the  binary  is  too  widely  separated  in  the  x-axis  for 
BINARY_FIT  and  because  the  secondary  is  positioned  beyond 
the  recorded  scan  in  the  y-axis.  In  the  NIRI  adaptive  optics 
image,  MT  59  has  a  companion  at  p  =  1"20  with  AK  —  2.75. 
This  total  separation  is  related  to  the  projected  separations  by 
=  Ax^  +  Ay^,  so  we  expect  that  secondary  would  have  a 
projected  separation  along  the  y-axis  of  Ay  =  — 0'.'91,  which  is 
beyond  the  scan  limits. 

MT  138.  We  determined  that  MT  138  (see  Figure  1.9)  is 
resolved  along  the  y-axis  with  Ay  =  -i-0'.'066  and  AV  = 
2.79  mag.  This  magnitude  difference  is  consistent  with  that 
from  a  close  companion  in  the  NIRI  image  at  p  =  1"5  and 
AK  =  3.18  mag.  According  to  the  FGS  aperture  position  angle 
(or  telescope  roll  angle),  the  NIRI  component  is  expected  to 
appear  at  Ay  =  -i-0'.'070  and  Ax  =  -i-l"338,  which  agrees  with 
the  BINARY_FIT  result  along  the  y-axis  and  puts  the  companion 
well  off  the  x-axis  scan. 

MT  304.  This  target  (Schulte  12  =  Cyg  OB2-12;  see  Fig¬ 
ure  1.18)  shows  evidence  of  a  very  close  companion  in  the 
CCF  and  second  derivative  tests,  with  values  of  a  =  8.81ct  and 
a  —  4.68ct  for  the  x-  and  y-axes,  respectively.  The  differential 
magnitude  from  the  BINARY_FIT  analysis  of  the  x-axis  scan 
was  used  as  a  constraint  in  the  fit  of  the  y-axis  scan,  to  arrive 
at  Ay  =  2.31  ±  0.21.  For  close  systems,  where  the  projected 
separation  is  less  than  the  size  of  the  fringe  (p  <  15  mas),  there 
is  an  ambiguity  in  the  “parity”  of  the  secondary  star’s  position, 
i.e.,  solutions  with  the  secondary  to  the  left  or  right  of  the  pri¬ 
mary  are  indistinguishable  from  one  another.  With  a  separation 
along  the  x-axis  of  61 .8  ±  3. 1  mas,  and  y-axis  solutions  yielding 
9.7  mas  or  —  14.9  mas,  the  position  angle  is  283.5  or  305. 9±  3.3 
degrees,  respectively.  We  were  not  able  to  detect  the  counter¬ 
part  in  the  NIRI  image,  because  the  separation  (p  =  63.6  mas) 
is  below  the  limiting  resolution  for  NIRI  (^^80  mas),  but  this 
companion  has  been  observed  using  other  interferometric  tech¬ 
niques  (R.  Millan-Gabet  2013,  private  communication).  MT  304 
is  an  early  B-type  hypergiant  of  very  high  luminosity,  but  it  may 
not  belong  to  the  class  of  Luminous  Blue  Variables  because 
its  flux  and  spectral  appearance  are  relatively  constant  (Clark 
et  al.  2012).  The  companion  we  And  is  too  faint  to  alter  the 
conclusion  about  the  star’s  high  luminosity.  If  we  assume  that 
the  projected  separation  corresponds  to  the  apastron  separation 
in  a  highly  elliptical  orbit,  then  the  orbital  period  is  P  30  yr, 
for  Ml  -I-  M2  =  120 Mq  (Clark  et  al.  2012)  and  d  —  \. A  kpc. 


Consequently,  additional  high  angular  resolution  observations 
over  the  next  few  decades  may  lead  to  a  mass  measurement  of 
this  extraordinary  star. 

MT  417.  This  star  (Schulte  22;  see  Figure  1.24)  is  the  only 
triple  system  resolved  in  our  sample,  and  we  developed  TRIPFIT 
to  model  the  scans.  Visual  inspection  of  the  y-axis  5-curve 
(central,  right  panel  of  Figure  1.24)  shows  that  there  is  a  very 
faint  third  component  near  Ay  =  -i-l"7.  Comparison  with  the 
NIRI  images  led  us  to  conclude  that  the  third  component  fringe 
is  blended  with  that  of  the  secondary  in  the  x-axis  scan.  The 
solution  for  differential  magnitudes  for  the  secondary  from 
the  two  axes  did  not  agree  with  each  other  (Amj^  —  0.85 
and  Arriy  —  0.45).  This  is  due  to  the  fact  that  the  second 
and  third  component  are  blended  in  the  x-axis  scan,  making 
the  amplitude  of  the  secondary’s  fringe  appear  smaller.  This 
is  not  the  case  for  the  y-axis,  where  the  two  components 
are  well  separated.  The  value  listed  in  Table  2  is  from  the 
y-axis,  which  has  the  larger  projected  separation,  and  the  error 
is  estimated  from  the  standard  deviation  between  the  fits  of 
the  different  calibrators.  MT  417  A  and  B  have  classifications 
of  03  If*  and  06  V((f)),  respectively  (Sota  et  al.  2011).  Our 
measurements  of  the  Ba,Bb  pair  (p  =  O'.'lOh,  9  —  160°,  and 
Am  —  2.48  mag)  agree  with  previous  works.  The  A,B  pair  was 
resolved  with  speckle  interferometry,  and  on  re-examining  the 
2007  speckle  data  using  the  methodology  of  Mason  et  al.  (2009), 
the  Ba,Bb  pair  were  observed  with  p  =  0'.'209,  9  =  18(T0.  This 
pair  was  also  resolved  through  AstraLux  and  HST  ACS/HRC 
imaging  by  Mafz  Apellaniz  (2010).  His  positions  and  magnitude 
differences  for  the  Ba,Bb  pair  (p  =  0'.'216,  9  —  181°48, 
and  Am  —  2.34  mag)  agree  with  the  speckle  data  and  our 
results,  within  the  uncertainties  due  to  the  blending  of  the  Ba,Bb 
components  along  one  axis  and  neglect  of  the  spatial  sensitivity 
of  the  photomultiplier  tubes. 

MT  429.  The  primary  of  MT  429  (see  Figure  1.25)  is  a 
short-period  eclipsing  binary,  and  the  presence  of  the  newly 
discovered  bright  third  companion  (AV  =  1.09  ±  0.02)  has  a 
strong  influence  on  the  interpretation  of  the  radial  velocity  and 
photometric  variations  (Kiminki  et  al.  2012). 

MT  516.  The  pair  of  MT  516  (see  Figure  1.38)  was  first 
resolved  through  speckle  interferometry  by  Mason  et  al.  (2009) 
who  designate  the  binary  as  WSI 67.  Our  position  and  magnitude 
differences  agree  within  the  uncertainties. 

MT  531.  MT  531  (see  Figure  1.39)  is  an  obvious  binary  in 
y  and  the  second  component  is  off  the  scan  in  x.  The  NIRI 
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image  shows  the  infrared  counterpart  with  p  —  1"45  and 
=  0.65  mag.  The  NIRI  position  predicts  that  the  companion 
has  projected  separation  along  the  jc-axis  Ax  —  — 1"38,  well 
beyond  the  recorded  scan. 

MT  605.  The  primary  of  MT  605  (see  Figure  1.45)  is  a 
double-lined  spectroscopic  binary  (Kiminki  et  al.  2012),  and 
the  unknown,  relatively  bright  third  companion  may  cause  line 
blending  difficulties  for  radial  velocity  measurements. 

MT  632.  This  system  was  an  obvious  binary  in  both  axes  (see 
Figure  1 .47).  We  resolved  the  secondary  in  the  NIRI  observation 
and  the  FGS  results  are  consistent  with  those  from  the  AO 
measurements. 

MT  696.  A  companion  to  this  star  (Schulte  27;  Figure  1.50) 
was  detected  using  the  second  derivative  test  (a  =  8.46a) 
and  the  cross-correlation  method  for  the  y-axis  scan,  (and  was 
detected  at  the  3a  level  in  the  second  derivative  test  of  the 
x-axis  scan).  We  adopted  the  differential  magnitude  from  the 
y-axis  solution,  because  of  the  greater  projected  separation  along 
the  y-axis.  There  exists  a  180°  ambiguity  in  the  position  of  the 
secondary  along  the  x-axis  (Ax  =  ±9.7  mas).  The  ambiguity 
is  reflected  in  the  two  solutions  for  the  position  angle  given  in 
Table  2.  The  primary  is  an  eclipsing,  double-lined  spectroscopic 
binary  (Kiminki  et  al.  2012),  and  the  presence  of  the  FGS 
companion  will  influence  the  interpretation  of  the  light  curve 
and  spectroscopy. 

SCHULTE  5.  This  target  (see  Figure  1.56)  is  an  obvious 
binary  in  the  scans  of  both  axes.  The  separation  is  too  large 
along  the  x-axis  for  a  BINARY_FIT  solution,  so  we  applied 
TRIPFIT  to  this  system.  The  program  arrived  at  similar  results 
for  the  differential  magnitude  derived  from  the  x  and  y  scans,  and 
we  list  the  average  in  Table  2.  The  uncertainties  were  determined 
from  the  standard  deviation  between  the  differential  magnitude 
from  both  axes.  This  system  was  first  resolved  by  Herbig  (1967). 
The  primary  is  a  short  period  {P  —  6.6  d)  eclipsing  binary 
consisting  of  two  luminous  evolved  stars  (Linder  et  al.  2009). 
Kennedy  et  al.  (2010)  suggest  that  there  is  another  star  close 
to  A  with  an  orbital  period  of  P  =  6.7  yr  based  upon  the 
variable  radio  emission.  Mason  et  al.  (2009)  did  not  resolve 
the  B  component,  probably  due  to  the  magnitude  difference, 
but  our  results  are  consistent  with  recent  observations  by  Mafz 
Apellaniz  (2010). 

6.  DISCUSSION 

Our  high  angular  resolution  survey  of  58  massive  stars  in  Cyg 
OB2  led  to  the  detection  of  13  resolved  systems  and  the  par¬ 
tial  resolution  of  two  other  stars.  The  resulting  binary  fraction 
of  22%  to  26%  is  consistent  with  the  results  from  Nelan  et  al. 
(2004),  who  resolved  5  out  of  23  OB  stars  (22%)  in  the  Carina 
Nebula  cluster.  We  were  able  to  find  the  infrared  counterpart  in 
our  NIRI  observations  for  1 1  of  the  resolved  systems,  but  not  for 
MT  304,  MT  696,  and  marginally  resolved  pairs.  Additional  ob¬ 
servations  with  ground-based  interferometry  confirm  the  com¬ 
panion  to  MT  304.  Two  of  the  most  massive  stars  in  the  associ¬ 
ation  (MT  304  =  Cyg  OB2-12  and  MT  417  =  Cyg  OB2-22A) 
have  resolved  companions  while  the  third  very  massive  star 
(MT  457  =  Cyg  OB2-7)  appears  single.  The  companion  of  the 
hypergiant  MT  304  (Cyg  OB2-12)  may  have  an  orbital  period 
of  a  few  decades,  and  continued  high  angular  resolution  obser¬ 
vations  should  reveal  the  companion’s  orbital  motion.  This  is 
potentially  a  very  important  target  for  mass  determination  (like 
HD  93129A;  Nelan  et  al.  2004). 

The  high  angular  resolution  capabilities  of  the  FGS  allow  us  to 
start  filling  in  the  observational  gap  in  the  period  distribution  of 


massive  binaries  (Mason  et  al.  1998).  The  angular  separations  of 
the  resolved  binaries  correspond  to  binary  orbital  periods  in  the 
range  of  20  <  P  <  20,000  yr.  The  longest  period  spectroscopic 
system  in  Cyg  OB2  observed  by  Kobulnicky  et  al.  (2012)  is  just 
over  6  yr,  only  a  factor  of  3  smaller  than  our  lower  limit.  This 
is  probably  the  smallest  gap  to  date  between  spectroscopic  and 
high  angular  resolution  methods.  Though  most  of  our  resolved 
companions  have  long  periods,  their  projected  separations  are 
less  than  10,000  AU.  This  corresponds  to  an  orbital  velocity 
that  is  larger  than  the  velocity  dispersion  of  the  association 
(a  =  8.03  km  s“';  Kiminki  et  al.  2008),  so  these  systems  are 
probably  orbitally  bound  companions. 

The  companions  we  detected  are  relatively  bright,  and  hence, 
it  is  important  to  account  for  their  flux  in  analyzing  the  spectra 
of  the  primary  stars.  Fainter  companions  are  undoubtedly 
present  in  this  period  range,  and  our  complementary  adaptive 
optics  study  with  NIRI  will  help  to  determine  the  mass  ratio 
distribution  of  lower  mass  and  fainter  companions  (at  least 
among  the  long  period  binaries).  The  systems  resolved  in 
our  sample  are  at  the  limits  of  what  is  angularly  resolvable 
with  a  single  aperture  telescope  today.  Preliminary  results 
from  the  larger  AO  sample  allows  us  to  perform  a  statistical 
study  of  chance  alignments,  which  suggests  that  the  systems 
resolved  with  the  FGS  have  <1%  probability  of  being  a  chance 
alignment.  In  a  future  work  we  will  combine  the  results  of  the 
spectroscopic  and  high  angular  resolution  surveys  of  Cyg  OB  2 
and  provide  an  unprecedented  census  of  the  binary  properties 
of  massive  stars  over  a  large  range  in  orbital  period. 
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